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ABSTRACT 

Context. The progenitor problem of Type la supernovae (SNe la) is still unsolved. Most of these events are thought to be explosions 
of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion scenarios, particularly those involving the externally 
triggered detonation of a sub-Chandrasekhar mass WD (sub-M C h WD), there is also a possibility of having an oxygen-neon (ONe) 
WD as progenitor. 

Aims. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with 
detonations in CO WDs of similar mass and observational data of SNe la. 

Methods. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses 
below the Chandrasekhar mass (M C h), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction net¬ 
work. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe la. We 
also perform binary evolution calculations to determine the number of SNe la involving ONe WDs relative to the number of other 
promising progenitor channels. 

Results. The ejecta structures of our simulated detonations in sub-M C h ONe WDs are similar to those from CO WDs. There are, 
however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are system¬ 
atically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO 
models. 

Conclusions. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should 
contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (~1.2M 0 ) shows that the less 
blueshifted spectral features fit the observations better, although they are too bright for normal SNe la. 

Key words, supemovae: general - nuclear reactions, nucleosynthesis, abundances - hydrodynamics, radiative transfer - white dwarfs 
- Stars: evolution 


1. Introduction 

Type la supemovae (SNe la) are believed to result from ther- 


ashes lead to loss of pressure support and collapse (Miyaji et al. 
1980; Miyaji & Nomotof 1987[i. They are therefore expected to 


monuclear explosions of white dwarf (WD) stars (Hoyle & 
Fowler 1960) in binary systems. The parameters of the progen 


form neutron stars rather than being disrupted in a thermonuclear 
explosion as their mass approaches Mch ( Nomoto| 1984| |1987[ 


Nomoto & Kondofl991jl. However, there is mounting evidence 


itor systems ( Wang & Han||2012), however, and the details of f or other progenitor channels contributing to (or dominating) 


the explosion mechanism (e.g. Hillebrandt & Niemeyer 


remain unclear. Several scenarios hold promise for exp 


normal or peculiar SNe la (see e.g. Hillebrandt et al. 


2000 ] ) 


aining 

2013] ). 


Typically, the exploding WD in all these scenarios is assumed 
to be a carbon-oxygen (CO) WD. Here, we explore the possibil¬ 
ity of explosions in ONe WDs. 

For a long time, explosions of near-Chandrasekhar mass 
WDs (near-Mch WDs) formed in the single-degenerate progeni¬ 
tor channel were the favoured model of SNe la (see [Hillebrandt 
& Niemeyer |2000 1 . In this scenario, ONe WDs are excluded 
because, although electron captures on 24 Mg and 20 Ne can ig¬ 
nite a nuclear flame, further electron captures in the O-burning 


the sample of SNe la (e.g. |Stritzinger et al. 

20061 |Ruiter et al. 

2009; Gilfanov & Bogdan 2010 Sim et al. 

2010] Scalzo et al. 

2014albl. Many of these alternatives involve detonations in sub- 


Mch WDs. This again raises the question of whether ONe WDs 
contribute to the progenitor population, since sub-Mch configu¬ 
rations are stable against gravitational collapse and detonations 
propagate rapidly enough such that electron captures do not lead 
to a collapse. This, however, requires the triggering of a deto¬ 
nation in ONe WD matter, which may be possible but has not 
yet been proven to work ( Shen & Bildsten| 2014i. Detonations 


in ONe WDs could be ignited in the double detonation (Fivne 
1990 Fivne & Glasner|[l990| 1991| Woosley & Weaver 1994 
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Finket al. 2007,2010; Woosley & Kasen(20TT, Moll & Woosley 

2013 

or violent merger (Pakmor et al. 2010 2012, 2013; Moll 

et al. 

2014) scenarios. Despite the uncertainties related to their 


ignition, we here investigate the question of how, if they do oc¬ 
cur, such events might differ from those with CO WD progen¬ 
itors and whether they might be identifiable as a subpopulation 
of SNe la. 


2. Population synthesis and the origin of ONe WDs 


Garcia-Berro et al.|( jl997| ) and |Gil-Pons & Garcia-Berro] ( |2001| ) 
investigated asymptotic giant branch (AGB) stars with zero-age 
main sequence (ZAMS) masses of 9-10 M G ; the first study is for 
a single star while the second paper describes the evolution of a 
close binary system. These stars produce ONe WDs with total 
masses of about 1.1 M Q . They are mainly composed of l6 0 and 
20 Ne, but they also contain some 12 C, the exact amount of which 
depends on the initial model Qpor a 9 M Q ZAMS progenitor, for 
instance, the l2 C mass fraction in the WD material can be up to 
~0.05, but it varies with radius. This C admixture is important 
for the initiation and subsequent propagation of the detonation 
because it acts as an accelerant for neon burning. However, suc¬ 
cessful initiation and propagation of the detonation are assump¬ 
tions in the present work. The viability of these can only be ad¬ 
dressed by spatially resolved direct numerical simulations of the 
hydrodynamics coupled to a full nuclear reaction network, un¬ 
fortunately still out of reach for full star explosion simulations. 

The WDs formed in binary systems cover a wider range 
of masses due to mass gain/loss from/to the binary companion 
during stellar evolution. Depending on previous mass transfer 
episodes, the ONe WD can easily have a mass ranging from 
1.08 Mo up to Mch. while CO WDs can be formed with masses 
as high as ~1.25M Q (see Hurley et al.] 2000|). We use the bi¬ 
nary population synthesis code StarTrack ( |Belczynski et ~ak 
2002, 2008| ) to predict the number of potential SN la progen 


itors that involve ONe WDs both for dynamical (mergers) and 
non-dynamical (classic double detonation) scenario^] and show 
their delay time distribution (DTD; Figure [T}. The results are 
summarized in Table [I] We find that averaged over a Hubble 
time, the total rate of CO+CO mergers is 1.06 x 10 13 M 0 1 
yr 1 , where the mass represents mass born in stars. This value is 
very close to the estimated SN la rate in Milky Way-like galaxies 
(1.1 x 10~ 13 M 0 -' yr -1 , see 


Badenes & Maoz 2012). Thus, 


we 

use the total number of all CO+CO mergers as a reference point 
when comparing total numbers for various progenitor scenarios 
in Table [I] 

If we compare the systems of ONe+ONe WDs combined 
with the CO+ONe WDs (denoted ONe+X mergers) that end up 
in a merger scenario, their fraction is 4% of the CO+CO mergers. 
Of course, not all CO+CO WD mergers involve primaries mas¬ 
sive enough to lead to a bright SN la in the violent merger sce¬ 
nario ( |Pakmor et al.|2010l|Ruiter et al.|2013] l. If we restrict the 
sample of CO+CO mergers to those where the primary WD has 
at least 0.9 M 0 , the ONe+X merger fraction is as high as 16%. 
We find that the number of binary systems involving ONe WDs 
that may lead to classical (non-dynamical mass transfer) double¬ 
detonations are about 8% of the classical CO double-detonations 


1 In addition to effects discussed here, another factor that influences 
the final composition of the AGB core is the carbon burning rate l[Chen| 


et al. 2014 1 . 

7 In the case of classic double detonations, the mass ratio is often 
sufficiently far from unity to enable stable mass transfer to proceed once 
the larger WD fills its Roche lobe. 


Table 1. Relative rates averaged over a Hubble time for theoreti¬ 
cally predicted explosion scenarios that may lead to SNe la. Data 
are extracted from the P-MDS model of |Ruiter et al.| l |2014) . All 
values have been normalized to the number of carbon-oxygen 
WD mergers (full mass range). The double-detonation systems 
are denoted by ddet. The accretion induced collapse systems 
arise from ONe WDs that collapse to a neutron star as they ap¬ 
proach Mch while accreting stably from a stellar companion. 
Entries marked with an asterisk are the systems for which we 
show delay times in Figure [T] 


SN la progenitor scenario rel. rate 

CO+CO mergers (all) 1.0 

CO+CO mergers (primary mass > 0.9 M 0 )* 0.27 

ONe+X mergers (all)* 0.04 

ddet (CO primary, all) 0.86 

ddet (CO primary mass > 0.9 M e )* 0.37 

ddet (ONe primary; all)* 0.03 

Chandrasekhar-mass CO WD (SD) 0.01 

accretion-induced collapse in ONe WD (AIC) 0.02 
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Fig. 1. Delay time distributions for a subset of SN la progenitors 
shown in Table [I] Top panel: mergers that involve one or more 
ONe WD (blue), and double detonations where the accretor is 
an ONe WD (orange). Bottom panel: mergers that involve two 
CO WDs where the primary (more massive) WD is > 0.9 M Q 
(blue), and double detonations where the accretor is a CO WD 
with mass > 0.9 M Q (orange). For both panels the combined 
DTD is shown in grey. 


(see Ruiter et al. 2014] ). Here, for CO WD primaries, we assume 
that only the systems in which the primaries are more massive 
than 0.9 M 0 will potentially lead to thermonuclear events that are 
bright enough to be considered SNe la (|Sim et al.|2010] |Ruiter| 


et al. 2014). The low-mass systems are thus excluded from being 


considered likely SN la candidates; however, we include their 
numbers in Table |T] for completeness. For the systems involv- 
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ing ONe primaries, the lowest mass primary is ~1.1M 0 , thus 
we include all of them. We note that for the population synthe¬ 
sis model presented, the number of ONe double-detonations is 
a factor of three higher than single degenerate (SD) systems in¬ 
volving a Mch CO WD. Taken together, these relative rate esti¬ 
mates suggest that potential progenitor systems in which an ONe 
WD explodes are frequent enough to constitute a substantial sub¬ 
set of thermonuclear explosions. 


In Figure [T] top panel, we show the DTDs for double deto¬ 
nations in ONe WDs and double WD mergers involving one or 
more ONe WD. The majority of ONe double detonations (or¬ 
ange) have delay times < 650 Myr. The donor stars in these bi¬ 
naries are naked helium-burning stars that were formed during 
the second common envelope phase that is encountered during 
the evolution. A small number of systems with delay times > 1 
Gyr involve helium WD donors. Most of the mergers (blue) are 
between ONe and CO WDs, though in some cases the mergers 
are ONe+ONe (14%). Taking the mergers alone, the DTD shape 
does not resemble a power law, which has been extensively ac¬ 
credited to merging CO WDs in the literature (see discussion in 
Totani et al.|2008]l. In the bottom panel of Figure[l]we show the 


DTD for binaries involving CO WDs where the primary CO WD 
mass for both double detonations and mergers is > 0.9 M Q . The 
double detonations show a very different DTD shape than that of 
the ONe systems; the peak at ~700 Myr is simply not found in 
systems with ONe primaries (instead there is a peak ~200 Myr). 
The lack of a later DTD peak for double detonations involving 
ONe WDs is due to the fact that events with delay times >800 


Myr typically involve degenerate (helium WD) donors (Ruiter 
et al. 2014) >, which are not produced as frequently in binaries in¬ 


volving the more massive ONe primaries. For mergers involving 
two CO WDs, again, the DTD shape is quite different in compar¬ 
ison to mergers involving an ONe WD. It is worth noting that the 
number of ultra-prompt (<100 Myr) mergers in binaries involv¬ 
ing an ONe WD is higher than that of CO+CO systems, even 
though the latter (even with the adopted mass cut) outnumber 
the ONe+X mergers by a factor of ~7. These ultra-prompt events 
originate from progenitors that undergo two common envelope 


events where the same star loses its envelope twice (see Ruiter 
et al. |2013| for discussion). Unlike the mergers with ONe WDs, 
the CO+CO merger DTD peaks around 200 - 300 Myr, consis¬ 
tent with the previous results of Ruiter et al. ( 2013| ). However, 
cutting out the lower mass CO WD primaries leads to a slightly 
flatter DTD shape than found when including a larger binary 
sample, where a t~ l power law is typically found. 


We do not include mergers between CO WDs and He-rich 
WDs, which may also lead to double detonations. We find that 
such systems always have CO WD masses below 0.7 M Q at the 
time of the merger. Typical masses for systems that lead to un¬ 
stable mass transfer and merge are <0.35 M 0 and ~0.6 M 0 for 
He and CO WDs, respectively. Whether or not a detonation of 
the CO core would be triggered by a He-detonation in such low- 
mass systems is still uncertain ( |Sim et al.|2012||Shen & Bildsten| 
2014[ >, and in fact some of these mergers may lead to the for¬ 
mation of RCrB stars rather than to thermonuclear explosions 
( |Webbink||1984| >. Either way, compared to SNe la, successful 
CO ignition in such systems would produce fainter and faster 
evolving thermonuclear transients (Sim et al.| |2012| l and would 
therefore not contribute to the SN la rate. 


3. Explosions of ONe white dwarfs 

3. 1. Numerical methods 

For our explosion simulations we use the Eulerian hydrody- 
namics code leafs (|Reinecke et al. 1999| 2002| Ropke 2005 [ 
Ropke & Hillebrandt 2 005 [i. This finite volume code is based 
on the Prometheus ( jFryxell et al. [1989 1 implementation of the 
“Piecewise Parabolic Method ’ (PPM, Colella & Woodward! 


1984). It includes an appropriate equation of state for WD matter 
based on the equation described by Timmes & Swesty (2000). 
The detonation front is modelled using the level set method 


(Reinecke et al.|| 1999t Golombek & Niemeyer |2005[ |Ropke & 

Niemeyer|2007| l, where the nuclear burning zone is numerically 
treated as infinitely thin. For numerical efficiency, instead of a 
full reaction network we use six pseudo-species, 12 C, 16 0, 20 Ne, 
4 He, intermediate mass elements (IME), and iron group ele¬ 
ments (IGE), approximately representing fuel and ash compo¬ 
sitions. Nuclear statistical equilibrium is followed by a tempera¬ 
ture and density dependent mixture of 4 He and the IGE species. 
The ash composition, which depends on fuel density and com¬ 
position, is read off from tables that are calibrated in a self- 
consistent manner as described below (Sect. [T2| >. 

This is sufficient to model the energetics of the detonation 
process. For convenience, the numerical simulations presented 
in this work are performed in two dimensions assuming axisym- 
metry. The numerical resolution is fixed to 1024 x 1024 grid 
cells that co-expand with the explosion of the WD so that all 
explosion ejecta can be followed to homologous expansion, ap¬ 
proximately reached 100 s after ignition ( Ropke|2005| l. To deter¬ 
mine the detailed chemical composition of the ejecta, we apply 
a nucleosynthetic postprocessing step (Travaglio et a l.||2004| >. It 
is based on approximately 90,000 tracer particles distributed in 
our 2D hydrodynamic al simulation. This is sufficient to obtain 
converged nucleosynthetic yields (Seitenzahl et al .|2010[ l. These 
tracer particles are passively advected with the hydrodynamical 
flow and record the thermodynamic trajectories of representa¬ 
tive fluid elements. In the postprocessing we follow the nuclear 
reactions, using a nuclear reaction network with 384 species 
( [Thielemann et al.| 1986|[19901 j 1996). 

To predict observables from our explosion models we con¬ 
duct radiative transfer simulations. Because of the spherical 
symmetry of the ejecta, we map the 2D distribution of the fi¬ 
nal chemical composition and the ejecta density to a ID grid in 
velocity space, using the same smooth-particle-hydrodynamics- 
like algorithm that is described in |Kromer et al. ( |2010|l. The ra- 
diative transfer calculations are performed with artis (|Sirnj2007 


Kromer & Sim||2009|l. We use the same atomic data as in |Gall| 


et al. 


(2012)) For the calculations 1.024 x 10 7 Monte Carlo pack¬ 


ets are used. The radiative transfer calculations start 2 days after 
explosion and end 120 days after explosion, with the simulation 
discretized into 111 logarithmically separated time steps. 


3.2. Calibration of the detonation model 

The detonation in our model is not resolved but is represented 
by a parameterized description. The energy release in the burn¬ 
ing process has to be determined and encoded in the fuel and 
ash composition with the six pseudo-species. To achieve con¬ 
sistent energetics, we employ an iterative calibration procedure 
( |Fink et al.|201Q[|Ohlmann et al.|2014] >. The tracer particles are 
arranged radially in the exploding star with constant distance in 
log p to ensure resolution of the transition regime between burn¬ 
ing to nuclear statistical equilibrium and incomplete Si burning. 
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Fig. 2. Illustration of the abundance table for the nuclear burn¬ 
ing. The composition of the ash, depending on the density for 
an initial composition of X( n C) = 0.03, 2f( 16 0) = 0.6, and 
X( 20 Ne) = 0.37. 

In an initial run, complete burning to IGE is assumed 
everywhere thus releasing the maximum amount of energy. 
Subsequently, a nucleosynthesis postprocessing step is per¬ 
formed that changes the ash composition in the lower density re¬ 
gions. The result is mapped into the six pseudo-species tabulated 
as a function of fuel density. This serves as input for the next it¬ 
eration step. The second parameter of our detonation model is 
the front propagation speed. It is set to the Chapman Jouguet 
velocity (CJ), which is determined according to the energy re¬ 
leased in each iteration step. For the calculation of the CJ speed, 
we use the equilibrium sound speed (Fickett & Davi s] 1979) . The 
procedure is repeated until the released energy in the explosion 
run matches the nucleosynthetic postprocessing result which is 
the case after ten iterations. The values of the calibrated table are 
plotted in Fig. [2] It clearly shows different burning stages, indi¬ 
cated by the pseudo-species. In the high density regime above 
~ 10 7 gcnT 3 fuel material is burned to nuclear statistical equi¬ 
librium (NSE) resulting in IGE after freeze-out. For intermediate 
fuel densities (10 6 < p [gem -3 ] < 10 7 ) the ash is composed of 
IME and oxygen. At lower densities, carbon and neon burn to 
oxygen, while below ~10 5 gem -3 burning ceases. 


3.3. Simulation set-up 


We calculate a series of explosion models of ONe WDs. These 
are set up in hydrostatic equilibrium with central densities po 
ranging (1... 2) x 10 8 g cm -3 , which is well below the threshold 


for ele ctron captures to become dynamically important dNomoto 


1987 


estimates p, 

8.5 x 10 9 g cm 


-3 


and 


Canal et al. 


1992 


9.5 x 10 9 g cm 

3 ). The initial temperature of all models is 


give 


assumed to be T = 5x10 5 K throughout the star. As in our previ¬ 
ous work on CO WDs (Sim et al. 2010) we assume uniform com¬ 
position with values motivated by the results of |Garcia-Berro| 
et al. (1997 1 and Gil-Pons & Garcla-Berro ( 2001] ). Specifically, 
the mass fractions of our initial composition are Jf( 12 C) = 0.03, 
X( le O) = 0.6 and Jf( 20 Ne) = 0.37. Our set-up procedure re¬ 
sults in ONe WDs with masses of 1.18 - 1.25 M G . Their param¬ 
eters are summarized in Table [2] Our simulations are for zero- 
metallicity main-sequence progenitors. Specifically, we do not 
include any intial abundance of the neutron rich isotope 22 Ne, 
which would slightly modify the results (Townsley et al.|2009| . 
As our intention is to study the outcome of detonations in ONe 
WDs rather than their progenitor evolution and ignition, we ig¬ 
nite the detonation by hand at the centre of the star. For compar¬ 
ison, we also run a simulation of a detonation in a CO WD set 



Fig. 3. 28 Si (upper panel) and 56 Ni (lower panel) mass fractions 
of the ejecta in velocity space. 

up with a central density of p = 1.5 x 10 8 g cm -3 , equal mass 
fractions of carbon and oxygen, and a uniform temperature of 
T = 5 x 10 5 K. 


3.4. Explosion simulations 

From our hydrodynamic al simulations we find many similari¬ 
ties, but also some clear differences between the CO WD and 
the ONe WD detonations. Table [2] shows the results in terms 
of kinetic energy of the ejecta and nuclear abundances of the 
pseudo-species. As expected from the differences in the binding 
energy of the fuel material relative to the ash (the energy release 
in burning 12 C to 56 Ni is about 3 0% higher than for burning 16 0 
or 20 Ne to 56 Ni; see Wang et al. 2012) the simulations lead to a 


significantly lower kinetic energy of the ejecta for ONe WD det¬ 
onations. The nucleosynthetic postprocessing results (Table [2]) 
show that there is a little less 56 Ni and overall IGE in the ejecta 
of the ONe WD detonation than in the ejecta of the CO WD 
explosion with the same initial mass. In contrast, the IME frac¬ 
tion of the ejecta is somewhat enhanced in the ONe detonation. 
Both effects are, however, not very pronounced. For the observ¬ 
ables, the more important difference is in the distribution of these 
species in velocity space. This is shown for 56 Ni and 28 Si as rep¬ 
resentative examples in Fig. [3] The distribution of both species 
is shifted towards lower velocities for the ONe WD explosion. 
Given that the 56 Ni masses are similar but the ejecta velocities 
are lower, we expect that the models will have similar brightness 
but the ONe WDs will give rise to smaller blueshifts of spectral 
features. Because of the lower expansion velocities, we expect 
an increase in the light curve rise time. 


3.5. Comparison of synthetic spectra and observations 

To compare our models with observations we have calculated 
synthetic spectra and light curves. Table [2] gives important val¬ 
ues from these calculations: the time of maximum brightness in 
B band the B-band decline time scale (Ani^(B)). and the 

magnitudes at fg max in the U , V. R, /, 7, II. K bands and in bolo- 
metric light. Fight curves for a subset of our models are shown 
(and compared to observations) in Fig. [4] There is relatively lit¬ 
tle variation among the light curves, as one would expect given 
that the masses of the models are similar. The difference in the 
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Table 2. Upper half: model parameters and most abundant nuclei of our nuclear network calculations. The first five models are 
for ONe WDs with a homogeneous composition of 3% 12 C, 60% 16 0, and 37% 20 Ne by mass. The last model, a homogenous 
CO WD with equal mass of 12 C and 16 0, is for comparison. Lower half: time of B-band maximum after explosion, absolute model 
magnitudes. Am 15 (B), and wavelength and corresponding blueshift velocity where the absorption of the Si 11 6355A feature is deepest 
at B-band maximum. 


Model 

ONelOe7 

ONel3e7 

ONel5e7 

ONel7e7 

ONe20e7 

C015e7 

Po [g cm-’] 

1.0 x 10 s 

1.3 x 10” 

1.5 x 10 s 

1.7 x 10 s 

2.0 x 10 s 

1.5 x 10 s 

Mo. [Mo] 

1.18 

1.21 

1.23 

1.24 

1.25 

1.23 

Bkin [10 M erg] 

1.14 

1.17 

1.19 

1.21 

1.22 

1.52 

IME [M 0 ] 

2.72e-l 

2.28e-l 

2.07e-l 

1.89e-l 

1.68e-l 

1,66e-l 

IGE [Mo] 

8.53e-l 

9.40e-l 

9.82e-l 

1.02e0 

1.06e0 

1.03e0 

lb O [Mo] 

5.20e-2 

4.06e-2 

3.52e-2 

3.12e-2 

2.65e-2 

2.24e-2 

24 Mg [Mo] 

7.21e-3 

5.72e-3 

4.98e-3 

4.41e-3 

3.76e-3 

6.10e-3 

28 Si [Mo] 

1.38e-l 

1.16e-l 

1.05e-l 

9.57e-2 

8.47e-2 

9.30e-2 

32 S [Mo] 

8.76e-2 

7.37e-2 

6.68e-2 

6.12e-2 

5.43e-2 

4.80e-2 

36 Ar [Mo] 

1.95e-2 

1.65e-2 

1.51e-2 

1.38e-2 

1.24e-2 

9.91e-3 

40 Ca [Mo] 

1.95e-2 

1.65e-2 

1.51e-2 

1.39e-2 

1.24e-2 

9.39e-3 

52 Fe [Mo] 

9.21e-3 

8.42e-3 

8.13e-3 

7.95e-3 

7.83e-3 

4.54e-3 

56 Ni [Mq] 

8.32e-l 

9.16e-l 

9.57e-l 

9.90e-l 

1.03e0 

l.OOeO 

57 Ni [Mo] 

6.43e-3 

7.72e-3 

8.39e-3 

8.98e-3 

9.74e-3 

1.19e-2 

r B m [days] 

19.1 

18.9 

18.7 

18.4 

18.4 

17.3 

m b 

-19.68 

-19.74 

-19.78 

-19.83 

-19.86 

-19.93 

Am 15(B) 

1.51 

1.47 

1.45 

1.40 

1.41 

1.67 

M U (t B m J 

-19.72 

-19.80 

-19.86 

-19.93 

-19.95 

-19.98 

M V^B m J 

-19.79 

-19.86 

-19.91 

-19.94 

-19.98 

-20.05 

M R (t B m J 

-19.57 

-19.62 

-19.64 

-19.65 

-19.67 

-19.69 

WbJ> 

-19.13 

-19.17 

-19.18 

-19.20 

-19.23 

-19.25 


-18.40 

-18.41 

-18.42 

-18.44 

-18.46 

-18.37 

M H<~ t B m J 

-17.76 

-17.80 
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-19.51 

-19.53 

-19.56 

Vsiii [ 10 y cm s '] 

1.3 

1.3 

1.4 

1.4 

1.4 

1.7 

'Li II [A] 

6090 

6080 

6060 

6060 

6060 

6000 


velocity structures of the CO WD and the ONe WD detonations 
lead to a rise time that is approximately two days longer for the 
latter. The time scales of the infrared light curve evolution are 
also slower in the ONe models (e.g. the secondary IR peak is 
reached a little later). 


Owing to their large 56 Ni masses, our models are too bright 
to provide a good match to a normal SN la such as SN 2005cf 

Fig# We do, how¬ 
ever, find better agreement with the light curves of SN 1991T 


(Pastorello et al. 2007; Garavini et al. 2007 


(Filippenko et al. 1992, Phillips et al.|1992 Ruiz-Lapuente et al. 


[1992;lira et al.: 1998). 

In Fig. [5] model spectra are shown at three epochs (-10, 0, 6 
days relative to fs m „). Overall, the spectra are similar and com¬ 
pare equally well to the observations; in particular, the impor¬ 
tant Si and S features are clearly present in the models with a 
strength comparable to the SN 2005cf data. There are, however, 
differences between the models in detail: 


- around maximum light the ONe WD detonation models 
show slightly stronger IME features; 

- the ONe WD explosion models have lower Si velocities 
(about ~ 400 km s" 1 ; see Table [ 2 J, compared to the CO WD 
explosion model. 

These results (see Table [2]i confirm what we expect from 
the hydrodynamical simulations (as discussed above). In all the 
models the blueshifts (e.g. of the Si 11 features) are generally too 
high compared to the observations, but this is less pronounced in 
the ONe WD detonations (Fig. [5]). 


The spectral features of our models are much too strong for 
SN 1991T. This holds in particular for the early epochs (-10 d, 0 
days after tB max ) where our models show strong Si 11 and Ca 11 fea¬ 
tures, while the spectra of SN 199IT do not show any strong ab¬ 
sorption lines. This discrepancy is consistent with expectations 
based on previous modelling of SN 1991T; in particular,|Sasdelli 


et al. ( |2014 1 have shown that the Si in SN 1991T should be pre 
dominantly located at velocities below ~ 12500 km s '. In con¬ 
trast, our models are Si-rich out to ~ 17000 km s 1 . Thus, despite 
having appropriate brightness, the spectra show that our models 
do not reproduce SN 1991T (or 91T-like objects in general). 


4. Conclusions 

Our population synthesis calculations show that the fraction 
of potential SN la progenitors involving ONe WDs is non- 
negligible. For the double-detonation scenario, they account for 
up to 3% of the sub-Mch systems (see Table [I]). For double de¬ 
generate mergers, those involving ONe WDs contribute 4%. If 
we restrict this census to systems that would reach sufficient 
brightness to produce a SN la event, this fraction increases to 
11%. In conclusion, explosions of ONe WDs are not a dominant 
channel of SNe la, but our results demonstrate that their con¬ 
tribution is important. It is therefore worthwhile to determine 
possible outcomes of thermonuclear explosions in ONe WDs in 
order to establish whether they can be identified as a subpopula¬ 
tion in transient surveys. 
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time (days since B lnax ) time (days since B max ) 


Fig. 4. Synthetic light curves in different bands for our ONe 


x 10 s g cm 3 and 2 x 


models with a central density of 1 

10 8 g cm~ 3 and our CO model compared to the observed 
light curves of SN 2005cf and SN 199IT. The SN_2005cf 
lightcurve is dereddened for E(B 


V ) = 0.097 ( |Past orello 
et al.||2007]l and the SN 199IT lightcurve is dereddened using 


E(B - V) = 0.13 ((Phillips et al.||1992|. The distance modules 
are 30.76 for SN 1991T (|Saha & Thmi|[2006| and 32.51 for 
SN 2006cf ( |Pastorello et al' 2007 1 . For the bolometric lightcurve 
of SN 199IT we used a NIR correction based on normal SNe la. 


X 

E 



4000 5000 6000 7000 8000 9000 


4[A] 


Fig. 5. Synthetic spectra for a subset of our detonation models as 
outlined by the labels in the bottom panel. For comparison, ob¬ 
served spectra of SN 2005cf (-10, 0, +6 d) and SN 1991T (-10, 
-3, +6 d) are shown. All times are relative to maximum bright¬ 
ness in B band. The spectra are de-redshifted: For SN2005cf by 
z = 0.00646 and for SN 1991T by z=0.006059. 

ing observed SNe la, the variation in the velocity distribution of 
species in the ejecta leads to changes in the predicted observ¬ 
ables, most importantly in the rise time of the light curves and 
Si line velocities (both are slower for ONe WD than for CO WD 
detonations). Thus, with respect to spectral line shifts, our ONe 
models fit the observations better than our CO models. In terms 
of brightness, however, all our models are too bright compared 
to normal SNe la. 

The question of whether particular explosion triggering 
mechanisms can be successful in ONe WDs remains open and 
has to be addressed in future work. 


We performed simulations of detonations in ONe WDs of 
different masses and compared these models to detonations in 
a relatively massive sub-Mch CO WD. As expected, the kinetic 
energy of the ejecta of ONe WD detonations is lower than that 
of an equal-mass CO WD. There is, however, very little differ¬ 
ence between the composition of the ejecta. Although overall 
there is no significant improvement or deterioration in match¬ 
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